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ABSTRACT 

We present the analysis of the X-ray brightness and temperature profiles for six clusters belonging to both 
the Cool Core and Non Cool Core classes, in terms of the Supermodel (SM) developed by Cavaliere, Lapi & 
Fusco-Femiano (2009). Based on the gravitational wells set by the dark matter halos, the SM straightforwardly 
expresses the equilibrium of the IntraCluster Plasma (ICP) modulated by the entropy deposited at the boundary 
by standing shocks from gravitational accretion, and injected at the center by outgoing blastwaves from merg- 
ers or from outbursts of Active Galactic Nuclei. The cluster set analyzed here highlights not only how simply 
the SM represents the main dichotomy Cool vs. Non Cool Core clusters in terms of a few ICP parameters gov- 
erning the radial entropy run, but also how accurately it fits even complex brightness and temperature profiles. 
For Cool Core clusters like A2199 and A2597, the SM with a low level of central entropy straightforwardly 
yields the characteristic peaked profile of the temperature marked by a decline toward the center, without re- 
quiring currently strong radiative cooling and high mass deposition rates. Non Cool Core clusters like A1656 
require instead a central entropy floor of a substantial level, and some like A2256 and even more A644 feature 
structured temperature profiles that also call for a definite floor extension; in such conditions the SM accurately 
fits the observations, and suggests that in these clusters the ICP has been just remolded by a merger event, in 
the way of a remnant cool core. The SM also predicts that dark matter halos with high concentration should 
correlate with flatter entropy profiles and steeper brightness in the outskirts; this is indeed the case with A1689, 
for which from X rays we find concentration values c 10, the hallmark of an early halo formation. Thus we 
show the SM to constitute a fast tool not only to provide wide libraries of accurate fits to X-ray temperature 
and density profiles, but also to retrieve from the ICP archives specific information concerning the physical 
histories of dark matter and baryons in the inner and the outer cluster regions. 

Subject headings: galaxies: clusters: general — galaxies: clusters: individual (A2199, A2597, A1689, A1656, 
A2256, A644) — X-rays: galaxies: clusters 



1. INTRODUCTION 

The evolution of a galaxy cluster comprises relaxed 
stretches punctuated by violent merger events affecting the 
gravitationally dominant dark matter (DM) halo, especially 
frequent in its early life after the initial collapse. 

A number of such events may be caught in action, as in the 
extreme case of the Bullet Cluster (Markevitch et al. 2002; 
Clowe et al. 2006). As to the many quieter clusters, one 
may wonder how long and how precise a memory of similar if 
less dramatic events — and of the more frequent ones caused 
by central outbursts from Active Galactic Nuclei (AGNs) — 
is retained by the hot diffuse baryons constituting their In- 
traCluster Plasma (ICP). We will show how a considerable 
amount of quantitative information may be retrieved even sev- 
eral Gyrs after the event from the X-ray bremsstrahlung emis- 
sion by the ICP in the inner cluster regions. 

On the other hand, the quiet stretches when external DM 
and gas are smoothly accreted across the virial radius R will 
shape the outer DM density and gravitational potential, and 
hence the ICP distribution and X-ray bremsstrahlung emis- 
sion in the cluster outskirts. Is there any simple and com- 
prehensive, physical modeling conducive to a precise yet fast 
analysis of both these regions? 

We will show how indeed our Supermodel (see Cavaliere, 
Lapi & Fusco-Femiano 2009, thereafter CLFF09) enables 
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US to accurately fit the full radial profiles of the X-ray ob- 
servables, namely, the temperature T{r) and the brightness 
S{r) oc n^(r) T^l^{r) that provides the baryon number density 
n. Whence we read out the imprints of the thermodynamic 
cluster history in the form of level, pattern, and time for de- 
positions or injections of ICP specific entropy k = k^T jn^^^ 
in Boltzmann units. Rather than relying on entropy profiles 
in the literature, we will focus on deriving the radial entropy 
runs from fitting the modeled primary observables for which 
precise, resolved and robust data (typically out to R/2) exist 
as is the case for the six clusters considered here. 

The Supermodel (SM) straightforwardly expresses the pro- 
files n{r) and T{r) in the entropy-modulated equilibrium of the 
ICP within the potential wells provided by the dominant DM. 
These two components are related not only by their common 
containing well, but also by parallel accretion of surrounding 
DM and baryons into the cluster volume. Thus the SM tells 
the thermal and dynamical past of a cluster back to its forma- 
tion time, to include subsequent mergers or AGN outbursts 
and current outer accretion. 

The SM is built upon the physical DM distributions derived 
by Lapi & Cavaliere (2009a) from the Jeans equation. These 
base on the radial run of the functional K = / p^l^, namely 
the 'DM entropy', that combines the 1-D velocity dispersion 
air) with the density p{r). From A^-body simulations Kir) is 
found to constitute a truly universal halo feature as it always 
follows a powerlaw run Kir) oc throughout the halos' bulk 
(Taylor & Navarro 2001; Dehnen & McLaughlin 2005; Hoff- 
man et al. 2007; Vass et al. 2009; Navarro et al. 2009). On the 
other hand, Lapi & Cavaliere (2009a) compute a, and find it 
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to lie in the range 1 .25 - 1 .3 from galaxies to massive clusters, 
narrowed down to 1.27- 1.3 from poor (overall mass M ~ a 
few 10^'^ Mq) to very rich (M ~ a few 10'^ M©) clusters. Two 
basic features of these 'a-profiles' entering the SM frame are 
briefly recalled in Appendix A. 

We stress that our treatment of the a-profiles (see Lapi & 
Cavaliere 2009a,b) agrees with recent numerical simulations 
(Zhao et al. 2003; Diemand et al. 2007; Wang & White 2008) 
in picturing the DM halo formation as a two-stage develop- 
ment. The epoch zt ranging from 0.5 to 2 marks the halo's 
transition from the stage of fast violent collapse scarred by 
major merging events, to one of progressively smoother and 
slower accretion building up the outskirts with little body per- 
turbation. 

During the latter, the DM profiles develop an increasing 
concentration c = R/r-2, that actually measures the current 
outer extension out to the virial radius R, relative to the re- 
gion inner to r_2 where p(r) is flatter than (see CLFF09). 
Starting with values c « 3.5 at the transition, c increases to 
current values c « 3.5 (1 +Zt), to attain values c « 10 or more 
for the 10% fraction of rich clusters with a transition as early 
as Z( « 1 .5 (see Lapi & Cavaliere 2009b); thus c measures the 
dynamic age zt of the cluster. Since R is typically of order 2 
Mpc, values c « 5- 10 imply r_2 ~ 100-200 kpc; in the fol- 
lowing, we will refer to 'outskirts' for the region r > r_2, and 
to 'inner body' for the region r < r_2. Based on the SM, we 
will discuss in particular how the outer distribution of the ICP 
entropy relates to the DM concentration. 

The SM describes the equilibrium of the ICP in the DM 
potential well, as we recall in § 2; there we briefly recap the 
main features of the ICP description in the form of the SM 
as developed by CLFF09, with the addition of Eq. (8) that 
enables computing the total mass. Then in § 3 we describe 
our analysis procedure based on the SM. In § 4 we apply such 
a procedure to the data for a set of six clusters with diverse 
profiles of X-ray brightness and temperature. Finally, in § 5 
we discuss the specific physical information we extract from 
our analysis of the X-ray observations. 

In our treatment we adopt a standard flat cosmology with 
normalized matter density SIm = 0.27, dark energy density 
J7a = 0.73 and Hubble constant = 70 km s"' Mpc"\ ex- 
cept when comparing with data otherwise scaled. 

2. THE SUPERMODEL 

We concentrate first on DM halos close to self-gravitational 
equilibrium under smooth and slow accretion, past the fast 
collapse stage and after any residual violent mergers, as dis- 
cussed in § 5. Within these DM gravitational wells, the ICP 
approaches hydrostatic conditions after outgoing blastwaves 
and shocks driven by central mergers or AGN outbursts have 
subsided. 

In such conditions, the ICP density is governed by the bal- 
ance between the gravitational force and the gradient of the 
pressure p = nksT/ij,; with the latter expressed as /? oc kn^^^ 
in terms of ICP entropy, the balance reads 



1 dikn^^^) GM{<r) 



fj,mpn 



dr 



(1) 



where /i « 0.6 is the mean molecular weight, ntp is the proton 
mass, and G is the gravitational constant. 

Once the radial entropy run k{r) is given as discussed be- 
low, the solution of this P' order differential equation allowed 
us (CLFF09) to write the profiles of the gas temperature and 



density in the form 



T{r)=k{r)ff-i^{r)=k^i\r) 



1 + 



dr' 



k-^l\f) 



(2) 

Here barred variables are normahzed to their boundary value 
at r =/?; the squared circular velocity v^(f) =M(< F)/F is taken 
from the a-profiles with their weak dependence on a, and is 
expressed in Appendix A in terms of hypergeometric func- 
tions; finally, we define bR = fimpV^(R)/kBT(R). 

The latter incorporates the boundary condition required for 
solving Eq. (1). It is physically convenient to fix such a refer- 
ence value at the virial radius r R rather than at the center, 
where v^(r) vanishes steeply (see Lapi & Cavaliere 2009a) 
while the ICP is often affected by violent stochastic events 
such as mergers and AGN outbursts. Moreover, at r ^ R 
closely universal accretion of external InterGalactic Medium 
(IGM) prevails for bofli Cool Core (CC) and Non Cool Core 
(NCC) clusters; this holds to lowest order, although the re- 
lated conversion of infall kinetic into thermal energy may dif- 
fer somewhat related to cluster age and preheating conditions, 
as discussed in § 3 and 5. 

In fact, in the absence of substantial preheating in the IGM 
larger than some 1/2 keV per particle (see Lapi et al. 2005; 
McCarthy et al. 2008), the energy conversion takes place in 
strong accretion shocks that linger at the virial radius (see 
Tozzi & Norman 2001; Lapi et al. 2005; Voit 2005), to imply 



(3) 



that takes on values around 2.5. These obtain from maximal 
conversion at a strong shock of the gravitational infall energy, 
that yields ksTR = |/imp A<I>r; this quantity involves the spe- 
cific kinetic energy Iv^AcpR gained by the IGM that starts 
from the turning point its free fall toward R across the poten- 
tial drop A^R = v\ A0R (see Lapi et al. 2005). The latter is 
provided by the DM a-profiles, and is conveniently normal- 
ized to the related v\ = v^iR); whence Eq. (3) follows. 

We stress that ksTR lowers when A(j> and the proportional 
infall kinetic energy are smaller owing to a large concentra- 
tion c; e.g., a value c « 10 holding for a cluster with an early 
transition (see § 1), in place of c « 4 holding for clusters with 
a recent transition, implies A0 to decrease from 0.57 to 0.47 
(see CLFF09). On the other hand, /cbTr also lowers when 
a high preheating level of the IGM weakens the shock and 
impairs the conversion efficiency as discussed by Lapi et al. 
(2005) and Voit (2005). 

To obtain the ICP temperature and density profiles the 
full radial run k(r) of the entropy is required. Its 
physical modeling is based upon the notions that en- 
tropy is eroded by radiative cooling on the timescale 
tc « 65ikBT/5 keV)'/2(n/10-3 crn^T^ Gyr (Sarazin 1988), 
whilst it is enhanced by shocks; at r « it is deposited by 
standing accretion shocks, and in the central region is injected 
by outbound blastwaves terminating in a shock, as are driven 
by supersonic outflows from AGN outbursts or head-on merg- 
ers. 

In the vicinity of r Ri /? the powerlaw approximation kocr" 
always applies (see CLFF09) with boundary value 



aR = 2.4-0.47 bR 



(4) 



varying around 1.1. Note that after Eq. (3) the upper bound 
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fls = 45/19 « 2.4 obtains for large as then ba tends to 
vanish; at the other end, decreases for weaker and weaker 
shocks corresponding to smaller Acf) and/or relatively stronger 
preheating levels. As no other major sources or sinks of en- 
tropy occur from the boundary at a few Mpcs down to the 
central r ^ 10^ kpc, in the outer range the entropy is deposited 
and stratified during the stage of slow accretion; thus a pow- 
erlaw radial run is set with slope a(r) that stays close to its 
boundary value ur. 

At the center, instead, entropy may be intermittently in- 
jected by shocks driven by mergers reaching down there (Mc- 
Carthy et al. 2007; Balogh et al. 2007), and by powerful 
AGNs residing in the central massive galaxies (see Valageas 
& Silk 1999; Wu et al. 2000; Scannapieco & Oh 2004; Lapi et 
al. 2005) as observed and reviewed by McNamara & Nulsen 
(2007) and Markevitch & Vikhlinin (2007). 

Thus the full entropy profile that combines central injec- 
tions with outer stratification may be described by the simple 
parametric expression (see Voit 2005 and references therein) 

kir) = kc + ii-kc)f^; (5) 

in fact, this approaches a constant value at small radii, and 
smoothly goes into a powerlaw at large radii. 

Entropy profiles similar to Eq. (5) have been recently re- 
ported by Cavagnolo et al. (2009) from an analysis of a Chan- 
dra archival sample comprising 239 clusters. In fact, most of 
them are well fit by a power law at large radii plus a constant 
value ^0 ^ kc at small radii, with a basically bimodal distribu- 
tion peaked at ko « 20 keV cm^ and at w 150 keV cm^. 

On the other hand, some clusters show evidence of a sharper 
entropy floor (e.g., Pratt et al. 2005 for A2218; see also Fig. 5 
in Cavagnolo et al. 2009) that we represent not only with 
a level kc but also with a definite extension r/, so that the 
corresponding radial entropy run reads 

k = kc (6) 

for F < Ff, and as 

k = k, + il-h)(^^^ (7) 

for r > Ff. The scale rj may be interpreted as the terminal 
radius just reached by an outbound blastwave driven by a vi- 
olent energy input at the center (see Lapi et al. 2005); at rj its 
decreasing Mach number Mir) has decayed to unity and the 
blast has stalled and degraded into adiabatic sound waves, as 
caught in action by Fabian et al. (2006) in the Perseus Cluster 
(see also § 5 for a discussion). In this picture, Eq. (5) repre- 
sents a later stage caused by diffusive smoothing and mixing 
of such a radial imprint, while the enhanced entropy level is 
stiU high before radiative erosion has set in. 

From the gas temperature and density profiles given by 
Eq. (2) we also derive the distribution of the gravitating mass 
in the form (see Sarazin 1988) 

M«P).3.65xlO- (^) (8) 

where the quantity rv^{r) grows slowly with r and saturates 
to unity. 

3. SUPERMODEL ANALYSIS 

Here we describe in detail how we use the SM to account 
for the X-ray brightness and temperature profiles for both 



classes of CC and NCC clusters (Molendi & Pizzolato 2001; 
Leccardi & Molendi 2009); in particular, we will examine 
in detail the following six clusters: A2199, A2597, A1689, 
A1656, A2256, A644. These have been collected from the 
literature on the basis of the quality and detail of the X-ray 
data, keeping a balance between the CC and NCC classes. We 
fit with the SM the profiles of emission weighted temperature 
and X-ray brightness, as given by Eqs. (9) and (10) below. 
The parameter values are pinned down on using a standard 
minimization procedure, and the uncertainties are quoted at 
the 68% confidence level. 

The free parameters are set as follows. The DM halo distri- 
butions depend weakly on the index a that we fix at the value 
a ~ 1.27, and are strongly marked by the concentration c that 
we leave as a free parameter. The ICP profiles are parame- 
terized by the slope a and the central value kc that define the 
radial entropy run after Eq. (5); in some clusters an accept- 
able fit requires to introduce the size rj- of the central floor 
after Eqs. (6) and (7). 

In the SM as with all models the virial radius R intervenes 
to set the data distance scale, as well as the bound to l.o.s. 
integrations; we keep it fixed when a robust determination 
is available from the literature, e.g., from observations of 
galaxy dynamics, or 'red sequence' termination, or gravita- 
tional lensing. Otherwise, we determine it from fitting (in the 
range where the data are reliable) the projected profile of the 
emission-weighted temperature 

f/^ d£ n\f) K[T{ry\ T{r) 
{ksTiw)) = kBTR _ , (9) 

dIn\r)Kmf)-\ 

in terms of the projected radius w = w/R. In the following 
we approximate the detailed cooUng function with A(r) oc 
T^l"^ as appropriate for the ICP in hot clusters. The above 
relation not only pins down the (horizontal) scale R, but also 
sets the normalization (vertical) scale Tr. The knowledge of 
R and Tr allows us to derive the ICP density Hr from fitting 
the brightness distribution 

x/ Un\r)^^l\r)F\E^,E^,^{J)■\■, (10) 

JO 

here 5o « 3.4 x 10"^^ (1 +z)~^ erg s~^ cm"^ arcmin"^, and the 
factor ^(£1 ,£2,7')- e"^'/*»^ -g-JhlksT ^^kes into account spe- 
cific instrumental bands E^-Ex (e.g., Ettori 2000). 
The SM actually predicts the values of Tr, ur and kR = 

kh^Rlnl^ from extrapolating the profiles into the outer clus- 
ter regions; observing the latter challenges the sensitivity and 
defies the resolution of most current instruments, but will con- 
stitute a main target for the next-generation X-ray telescopes 
planned to study low surface brightness plasmas (see § 5). 

With the use of such facilities R will be read out from the 
profiles as the position of the shock discontinuities in n{r) and 
T{r) above the values prevailing in the IGM (see Lapi et al. 
2005), even though such discontinuities may be blended by 
complex texture of the shocks and smoothed by projection 
integrations (see Tormen et al. 2004). On the other hand, 
substantial preheating of the IGM may weaken the accretion 
shocks and lower Tr. Such conditions will be indicated by any 
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Fig. 1 . — Abell 2199. Left panel: The solid line is our fit (x^ = 73.9/56) to the brightness profile measured by Mohr et al. (1999), on adopting the entropy 
profile given by Eq. (5). Right panel: The solid line is our fit (x^ = 15.8/5) to the radial temperature profile measured by Snowden et al. (2008), on adopting the 
entropy profile given by Eq. (5). 



discrepancy between the value resulting from Eq. (9) and the 
reference value ksTR corresponding to strong shocks in a halo 
of given concentration. 

4. RESULTS FOR INDIVIDUAL CLUSTERS 

Next we present the results for each of our six clusters; the 
individual fitting parameters are given in the corresponding 
subsections and are collected in Table 1. 

4.1. Abell 2199 

A2199 is a rich cluster at z « 0.031 centered on the cD 
galaxy NGC6166 that hosts the radio source 3C338. The lat- 
ter shows evidence of several events of radio emission, as wit- 
nessed by the clear correspondence between positions of radio 
lobes and depressed X-ray surface brightness ('cavities'). 

This is a CC cluster with the temperature profile declin- 
ing from about 200 kpc toward the center as observed with 
Chandra by Johnstone et al. (2002) and with XMM-JVewton 
by Snowden et al. (2008). Fig. 1 shows the SM fits to the 
projected emission- weighted temperature profile given by the 
latter authors, and to the brightness distribution by Mohr et 
al. (1999) obtained with ROSAT. Both profiles do not show 
evidence of an entropy floor down to r « 2 kpc. 

From the SM, we obtain the virial radius R = 2.1;|^4 Mpc 
and the halo concentration parameter c = 6. ?;!;[{]. In the inner 
ICP regions we find h = 0.39^^ }^ x lO'l Throughout the 
cluster body we derive the entropy slope a = 0.95;|^ qJ. In the 
outskirts we find keTR = 1-92,^% keV and hr = L40!g:g^ x 
10"^ cm"^, yielding kR = 3320;!;}3o keV cm^. Correspondingly, 
we obtain kc = 13^g keV cm-^. 

4.2. Abell 2597 

A2597 is a nearby cluster at z w 0.08 that harbors the cen- 
tral radio source PKS 2322-123; it shows features in dif- 
ferent spectral regions. In particular, Chandra X-ray obser- 
vations show a bright central region and two low surface 
brightness ghost cavities (McNamara et al. 2001), adding to 
other core structures suggestive of irmer cavities (Clarke et 
al. 2005). A2597 is a CC cluster as is also confirmed with 



XMM-Newton observations by Morris & Fabian (2005) and 
Snowden et al. (2008). The average temperature is around 2.6 
keV (Sarazin & McNamara 1997). 

We fit (see Fig. 2) the surface brightness data obtained 
by ROSAT/PSPC observations (Xue & Wu 2000) and the 
temperature profile observed by XMM-JVewton (Snowden et 
al. 2008), on considering the two entropy profiles given by 
Eqs. (5) or by Eqs. (6) and (7). The latter include the addi- 
tional parameter r/ = 50^7 kpc that turns out to be significant 
at the 99.97% level after the F-test, and may be interpreted as 
the effect of the central radio source PKS2322-123. On the 
other hand, the size of the entropy floor lies within the first 
bin of the temperature data (see Fig. 3), below the resolution 
of our profile. 

From the SM, we obtain the virial radius R = i-9'^)'l Mpc 
and the halo concentration parameter c = 7.2;|;5 2. In the inner 
ICP regions we find h = 0.21^g:|| x 10"^. Throughout the 
cluster body we derive the entropy slope a = 0.71^Q{5g. In the 
outskirts we obtain kgTR = 2. l!g;^ keV and hr = 2.0!j] ;| x 10"^ 
cm~^, yielding kR = 2850:1:94™ keV cm^. Correspondingly, we 
obtain kc = 6tl^ keV cm^. 

4.3. Abell 1689 

A1689 is a rich cluster at z « 0.183 that has been studied 
in various spectral regions; a joint analysis of HST/ACS and 
Chandra measurements has been carried out by Lemze et al. 
(2008) to derive the profiles of the ICP density and tempera- 
ture. The cluster is centered on a cD galaxy, is roughly spheri- 
cal and appears fairly relaxed (but see Andersson & Madejski 
2004). The X-ray luminosity is Lx « lO"'^ erg s"^ in the en- 
ergy band 0.5-7 keV, due to the hot ICP with kgT « 9.4 keV 
(Ebeling et al. 1996; Xue & Wu 2002; Andersson & Made- 
jski 2004). Strong and weak lensing observations (see Lemze 
et al. 2008 and references therein) show the projected mass 
profile to continuously flatten toward the center, with a steep 
outer profile that require a high concentration c w 13.7:l:[ | (see 
Broadhurst et al. 2008 and references therein), the hallmark 
of an early transition. 

With the SM we have fitted the ICP temperature profile ob- 
served with XMM-Newton by Snowden et al. (2008), and the 
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Fig. 3. — A2597. The solid line is our fit (x^ = 1.1/4) to the temperature 
profile measured by Snowden et al. (2008), on adopting the entropy profile 
given by Eq. (5). 

surface brightness profile observed by Mohr et al. (1999) and 
by Lemze et al. (2008). In agreement with the latter authors, 
our fits (see Fig. 4) do not require an entropy floor; any central 
flattening may affect only radii r < 2 x 10~^R « 4 kpc. 

We adopt the virial radius R = 2.l Mpc fixed at the value 
given by Lemze et al. (2008), and from the SM we find the 
halo concentration parameter c = 13.6^1. In the inner ICP 
regions we find kc « 2.4;^^ g x 10"^. Throughout the cluster 
body we derive the low empirical value of the entropy slope 
a = 0.80:!:[] "6 (consistent with the value 0.82^° measured by 
Lemze et al. 2008), that we relate after Eqs. (3) and (4) to the 
high concentration. In the outskirts we obtain ksTR = 4.4;!^ g 
keV and = 3.2!|; ;[ x lO'^ cm-^ yielding yt« « 4360!gg keV 
cm^. Correspondingly, we obtain kc = I05'!i'll cm^. 

Finally, in Fig. 5 we show the mass profile derived from the 
SM after Eq. (8) on using the parameters obtained from our 
fits; this turns out to be in good agreement with the total mass 



profile obtained from the lensing data by Lemze et al. (2008). 
4.4. Abell 1656 (Coma Cluster) 

A 1656, the popular Coma Cluster, is a very rich cluster at 
Z « 0.023. Its ICP at temperature kgT « 8 keV is a powerful 
X-ray emitter with luminosity Lx ^ 10"*^ erg s~' (see David 
et al. 1993) extending out to about 1° from the center (e.g., 
Briel et al. 2001). 

Several determinations of the virial radius ranging between 
2 and 3 Mpc are given in the literature (Castander et al. 2001 ; 
Lokas & Mamon 2003; Kubo et al. 2007; Gavazzi et al. 
2009). For our analysis, we have adopted the value of 2.2 
Mpc reported by Gavazzi et al. (2009), and checked that our 
results depend only weakly on this choice within 1 standard 
deviation. 

From the SM we find the halo concentration parameter 
c = 3.0;|;Q g, the hallmark of a recent formation (see § 1). In the 
inner ICP regions we find kc « 10^ j x 10"^. Throughout the 
cluster body we derive the entropy slope a = 1.30;!;o 23. In the 
outskirts we obtain ksTR = 5.7+} j] keV and Hr = 3.66!^;}^ x 
10"^ cm-^ yielding kR = 5170![!]^^ keV cm^ Correspond- 
ingly, we find kc = 520![^j] keV cml 

Our fits to the the brightness data observed with ROSAT 
by Mohr et al. (1999) are obtained on using the two entropy 
profiles given by Eqs. (5) or by Eqs. (6) and (7), and are illus- 
trated in Fig. 6; the values in the caption strongly indicate 
the presence of an entropy floor with extension r/ = 250;!^4 
kpc. On the other hand, the fit to the emission- weighted tem- 
perature profile (see Fig. 7) obtained from XMM-Newton ob- 
servations by Snowden et al. (2008) is roughly isothermal, 
and within its resolution does not require by itself an entropy 
floor 

All the above concurs to a picture of A 1656 as a halo that 
has just collapsed and undergone major mergers, as confirmed 
by the structured features around the main galaxies NGC4874 
and NGC4889 and by the ongoing faU of the NGC4839 
galaxy group into the main body of the cluster (Adami et 
al. 2005). A1656 also exhibits extended radio emissions, in- 
cluding a young radio halo close to center and an outer radio 
relic located in the SW direction beyond the NGC4839 galaxy 
group (Giovannini et al. 1991, and references therein). 
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Fig. 4. — A1689. Left panel: The solid line is our fit (x^ = 19.8/49) to the brightness profile measured by Lemze et al. (2008), on adopting the entropy profile 
given by Eq. (5). Right panel: The solid Une is our fit (x^ = 9.2/2) to the temperature profile measured by Snowden et al. (2008), on adopting the entropy profile 
given by Eq. (5). 
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Fig. 5. — A1689. The soUd line is our fit to the mass profile measured by 
Lemze et al. (2008) through gravitational lensing observations. 

From Eq. (8) we deiive an overall mass Mr = 1.24'^-^ x 
IO'^Mq, in agreement with the value of 9.7+3 ' x IO'^'M© ob- 
tained by Gavazzi et al. (2009) from gravitational lensing. 

4.5. Abell2256 

A2256 at z « 0.06 is a complex cluster, as highlighted by 
several studies in various spectral bands; it is a strong X- 
ray emitter with luminosity Lx « erg s"'. We have used 
the SM to fit the ROSAT brightness distribution by Mohr et 
al. (1999) and the temperature profile observed with XMM- 
Newton by Snowden et al. (2008), on using the entropy pro- 
files given by Eqs. (5) or by Eqs. (6) and (7); our results are 
illustrated in Fig. 8 and 9. 

From the SM, we find the virial radius R = 2.2+° 3 Mpc, and 
the halo concentration parameter c = 2.7"""'^. In the inner ICP 
regions we find kc = 6.2+3 j x 10"^. Throughout the cluster 
body we derive the entropy slope a=l .48!^ 29. In the outskirts 



we obtain keTR = 4.4![;| keV and hr = 3.64^^^ x lO'^ cm'^, 
yielding kR = 4000^1^1'' keV cml Correspondingly, we find 
kc = 24^tjlt keV cm^. 

Both the brightness and temperature distributions call for 
an entropy floor; basing on Eqs. (6) and (7) we derive a floor 
radius rj = 264+Jq^ kpc from the brightness, an evidence re- 
inforced by the value r/ = 265+J7Q kpc we obtain from the 
temperature profile. Its introduction allows the SM to fit weU 
the structured temperature profile of A2256; this features a 
temperature decrement similar to a CC cluster (e.g., Piffaretti 
et al. 2005; Leccardi & Molendi 2009), but at small radii T(r) 
reverses its trend and increases toward the center. Such a be- 
havior is understood from the relation T(r) oc k{r)n{r)^l^; in 
the inner region where the entropy is constant, the tempera- 
ture is expected to decrease outwards foUowing the density, 
while for r > r/ the entropy starts to increase and to domi- 
nate the density decrement, so raising the temperature out to 
a peak at r « 350 kpc (about 5'); beyond the peak the density 
steepens and offsets the entropy rise. The central temperature 
behavior suggests that the energy delivered by a merger has 
remolded the whole inner structure, and hence that the ICP is 
itself thermodynamically young within a dynamically young 
DM halo. 

The halo's young age is supported by the low values c « 4 
of the DM concentration, as expected in clusters with a re- 
cent transition from fast coUapse to slow accretion (see § 1). 
It is also indicated by the recent and intense merger activity 
that characterize A2256, as pinned down by ROSAT, Chan- 
dra and XMM-Newton observations (Briel et al. 1991; Sun 
et al. 2002; Berrington et al. 2002; Miller et al. 2003; Bour- 
din & Mazzotta 2008). Such a merger activity is quite recent 
as indicated by the presence of relativistic electrons with life- 
time of 10"' Gyr that enlighten the radio halo and the very 
extended, bright relic in the NW region of the cluster (Bridle 
& Formalont 1976; Kim 1999; Clarke & Ensslin 2006). All 
that concurs with the high central entropy level pinned down 
by the SM analysis. 



4.6. Abell644 
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Fig. 6. — A1656. Left panel: The solid line is our fit (x^ = 408.5/44) to the brightness profile measured by Mohr et al. (1999), on adopting the entropy profile 
given by Eq. (5). Right panel: The solid line is our fit (x^ = 53.0/43) on adopting the entropy profile given by Eqs. (6) and (7). 
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FIG. 7. — A1656. Left panel: The soUd line is our fit (x^ = 3.7/6) to the temperature profile measured by Snowden et al. (2008), on adopting the entropy profile 
given by Eq. (5). Right panel: The solid line is our fit (x^ = 3.7/5) on adopting the entropy profile given by Eqs. (6) and (7). 



A644 is a radio quiet cluster at z w 0.07. The X-ray imag- 
ing with Chandra by Buote et al. (2005) shows that, while the 
outer ICP is smooth and relaxed, the central ^ 10^ kpc fea- 
tures very complex and interesting conditions. Although the 
latter do not bear strong marks of past radio activity like X-ray 
cavities or filaments, the peak of the X-ray emission is found 
to be shifted by about 60 kpc from the cluster centroid (see 
also Bauer & Sarazin 2000), beyond the cD galaxy located at 
about 40 kpc. The temperature profiles differ when centered 
on the centroid or on the X-ray peak; the former profiles de- 
viates only weakly from a constant, while the latter shows a 
behavior similar to, if sharper than A2256, in that it decreases 
considerably toward the center before reversing its course and 
rising again at small radii (Buote et al. 2005). Our SM fits to 
this temperature profile, obtained on using Eq. (5) or Eqs. (6) 
and (7), are illustrated in Fig. 10; the fit to the brightness pro- 
file centered on the cluster centroid is illustrated in Fig. 11. 

We find the virial radius R = 2.l'^'l Mpc and the halo 
concentration parameter c = 3.9-0.2. At the center we find 



kc = O.lt^ l X 10 ^. On using Eqs. (6) and (7) we derive a ra- 
dius rf = 61^4^1 kpc, whose need is substantiated by the F-test 
at 99.9% level and also strengthened by the fit to the Chandra 
brightness profile centered on the X-ray peak (D. Buote, pri- 
vate communication) that consistently yields rf = 66tl kpc. 
Throughout the cluster body we derive the entropy slope 



a 



1.06_|J jj. In the outskirts we find ksTR 
Hr = 2M^ \l X 10-5 cm-\ yielding kg = 5200!^}^ keV cm^ 



4.9^{| keV and 



36;|;?A keV cm^. Finally, from 

■)+0.7 



""-0.19 ^ 

Correspondingly, we find kc - -'^-20 
Eq. (8) we derive the total mass M = 1.2+[J;7 x lO'^/W© 

Note that if one insisted on applying the SM also to temper- 
ature and brightness profiles from the centroid on the basis of 
Eqs. (6) and (7), one would obtain r/ = 104;!^ kpc and a related 
lower bound kc w 124^'^" keV cm^. The large variance in the 
above parameters entering the entropy floors signals complex 
substructures, that may be interpreted as a high density, low 
entropy clump ('cold drop') around the X-ray peak. Our re- 
sults agrees with the analysis by Buote et al. (2005), who 
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Fig. 8. — A2256. Left panel: The solid line is our fit (x^ = 434.1/65) to the brightness distribution measured by Mohr et al. (1999), on adopting the entropy 
profile given by Eq. (5). Right panel: The solid line is our fit (x^ = 109.7/64) on adopting the entropy profile given by Eqs. (6) and (7). 




r (arcmin) 



r (arcmin) 



Fig. 9. — A2256. Left panel: The solid line is our fit (x^ = 12.4/5) to the temperature profile measured by Snowden et al. (2008), on adopting the entropy 
profile given by Eq. (5). Right panel: The solid line is our fit (x^ = 3.6/4) on adopting the entropy profile given by Eqs. (6) and (7). Note that the first bin covers 
the range out to about 0.05 /ij^ Mpc in Fig. 8. 



in tenns of two differently centered /3-models (Cavaliere & 
Fusco-Femiano 1976) find two core sizes similar to our ex- 
tensions rf. 

In view of the lack of radio emission and X-ray cavities, 
this complexity may be understood in terms of a merger hav- 
ing just remolded an iimer region of the ICP (see Henning et 
al. 2009). We stress that such ICP substructures constitute a 
common, but progressively more pronounced trait, of A1656, 
A2256 and A644. 

5. DISCUSSION AND CONCLUSIONS 

In this paper we have analyzed with the Supermodel (SM) 
the profiles of X-ray temperature and surface brightness of 
the IntraCluster Plasma (ICP) in a set of six clusters (adding 
to the three ones preliminarily reported in CLFF09) with ex- 
isting detailed data. We have shown how effective is our SM 
to represent and understand the main Cool Core /Non Cool 
Core dichotomy in terms of two physical parameters marking 
the full ICP entropy profile: the central value kc, and the outer 



slope a (see Figs. 1-7). Moreover, the SM makes sense of 
more structured profiles (see Figs. 8-11) in terms of the ad- 
ditional, physical parameter rf marking the extension of the 
entropy floor. 

The working of the SM may be reduced to the bones as 
follows. The spatial scale for the temperature peak in the 
CCs, and for the outward temperature and density declines in 
all clusters, is set by the underlying gravitationally dominant 
DM distribution; specifically, such a scale is set by the peak 
of the DM velocity dispersion a^{r) at the position « r_2, 
that divides the 'inner' from the 'outer' cluster regions (see 
Figs. 1 and 3 of CLFF09). In fact, the SM ensures that the ap- 
proximation T ^ (T^/Aii [with /?„, = i.impCp'{r„)/kBT{rm), see 
CLFF09J is to hold closely for all clusters around and over 
a wide radial range for the CCs"*. In the inner ICP regions, the 
temperature and density profiles are governed primarily by the 

* The formal reason why around the temperature r(r) mirrors so well 
the DM dispersion <T^(r) is that the latter can be expressed quite similarly to 
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Fig. 10. — A644. Left panel: The solid line is our fit (x^ = 15.2/5) to the temperature profile centered on the X-ray peak measured by Buote et al. (2005), on 
adopting the entropy profile given by Eq. (5). Right panel: The solid line is our fit (x^ = 1.97/4) on adopting the entropy profile given by Eqs. (6) and (7). Note 
that the first bin covers the range out to about 0.04 /ij^ Mpc in Fig. 11. 




Fig. 1 1 . — A644. Left panel: The solid line is our fit (x^ = 120.7/56) to the brightness profile measured by Mohr et al. (1999), on adopting the entropy profile 
given by Eq. (5). Right panel: The solid line is our fit (x^ = 76.3/55) on adopting the entropy profile given by Eq. (6) and (7). 



central entropy level kc set by thermodynamical events, i.e., 
energy discharged and blasts driven by major mergers or AGN 
outbursts; in a number of cases the latter imprint substructure 
on the small scale r/ comparable to r„,. A novel feature of 
the SM (relative to handy isothermal or polytropic /3-models 
where the entropy is assumed to be a functional k oc n^~^^^ of 
the density, see Cavaliere & Fusco-Femiano 1976) is consti- 
tuted by the articulated radial entropy run entering Eq. (2). 

We have collected in Table 1 the SM fitting parameters for 
all clusters in our set. By inspection it is apparent a correlation 
between the inner ICP profile type (marked by the CC/NCC 
tags) with the central entropy level kc, the outer entropy slope 
fl, and the DM concentration c; high values of c and low val- 
ues of a and correspond to the CC class, while the opposite 
trend holds for NCC. We understand these trends in the frame- 
work of two-stage cluster formation (see § 1) as follows. For 

Eq. (2), except for the constant term within square brackets that is anyway 
negligible in cluster bodies (see CLI'F09), and particularly so for the CCs. 



example, low values of a correspond to high values of Br (see 
Eq. 4) owing to low values of (see Eq. 3); these are related 
to high concentrations c= 3.5(1 +Z() (see § 1 and 2), which 
imply early transition redshifts z,, i.e., an old age. We stress 
that in the six clusters considered here the outer parameters 
a and c turn out to be related to the inner parameter kc (and 
rf when applicable) in terms of the cluster age, as explained 
below. 

Thus the CC clusters with their low value of a appear to 
be generally older structures, currently in their stage of slow 
and smooth accretion, with the main action taking place in the 
outskirts under the form of calm entropy deposition by gravi- 
tational accretion shocks. Toward the center, the CC hallmark 
is constituted by a temperature peak overlapping (cf. Fig. 3 
in CLFF09) the peak in at r,„ « 10"'/? (cf. the profiles of 
A2199 and A2597 in Figs. 1-3). The condition for the peak 
to occur after the SM is a low value of the central entropy 
kc<3x 10"2, which comes to 30-50 keV cm^. This behav- 
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TABLE 1 

Fitting parameters from the SM analysis 



Cluster 


Class 


c 


a 


kc [10-2] 


r/ [kpc] 


R [Mpc] 


^ijTii [keV] 


A2199 
A2597 
A1689 
A1656 
A2256 
A644 


CC 
CC 

cc 

NCC 

NCC (RCC) 
NCC (RCC) 


6.7+'° 
3.9-0.2 


95+O 01 
-Q.Ql 

' '-0.05 

J 3q+B;5H 


39+0.16 

21+^-^^ 
2 4+0.8 

10+! 

6.2+3' 
7+^:1 

"•'-0.1 


<2 

50+? 

<4 

25^!^ 
254+102 
^ -80 

66!| 


2.1-0.4 

1 Q+0.4 
^•^-0.4 

2.1 

2.2 

2 2+*-3 
2 1+8-^ 


1 03+0.05 

4 4+8:^ 

5.7+^'^ 
44+i:8 

4 Q+?:^ 



Note. — The classification is taken from Molendi & Pizzolato (2001), see also Henning et al. (2009). For A1689 and A1656, the value of the virial radius is 
taken from the literature (see text). The values of ksTR are computed from Eq. (9), and are to be compared with those from the strong shock condition (see § 2 
and 3). For A644 the values of kc and r/ refer to the X-ray peak region; this explains the low kc value (see Fig. 10). 



ior is highlighted in terms of T{r) oc k{r)n{rf-/^; as the ICP 
density n{r) rises monotonically inward, T{r) will peak and 
then decline toward the center as k{r) decreases sharply to- 
ward a low central value kc. 

Such an inward decline of T( f) to a low but finite central 
value Tc oc A;" "'^ links to a high density oc to consti- 
tute the cool core, a feature of the non-radiative SM equilib- 
rium (also present in simulations discussed by Borgani et al. 
2008). As expanded upon by CLFF09, the SM does not in- 
clude enhanced cooling, even less any related inflow; it rather 
focuses the conditions for enhanced radiation and fast cool- 
ing to set in on the timescale tc « 0.3(A:c/15 keV cm^)'-^ Gyr. 
This would lead to a cooling catastrophe (e.g.. White & Rees 
1978; Blanchard et al. 1992), that may be stabilized by ICP 
condensing around and into a central massive galaxy to trigger 
accretion on the nuclear black hole. These conditions kindle 
up AGN activities that drive rising bubbles or even outgoing 
blastwaves, feed back entropy, and distribute it widely into 
the ICP (see Binney & Tabor 1995; Ciotti & Ostriker 2001; 
CavaUere et al. 2002; Churazov et al. 2005; Lapi et al. 2005; 
Voit & Donahue 2005; Tucker et al. 2007). 

We have analyzed in detail the two CC clusters A2199 and 
A2597 with their inward decrease of the temperature. We 
have found low central entropy levels kc'^\5 keV cm^ typical 
of CCs, with little or no need for an extended entropy floor. 
We have derived outer powerlaw slopes a < 1 (see Table 1), 
lower than the standard value 1 . 1 corresponding to the stan- 
dard concentration c « 4. 

This trend culminates with A 1689, a cluster with a CC-like 
inner profile but featuring interesting outer pecuharities. On 
an empirical stand, our SM analysis confirms the results by 
Lemze et al. (2008) concerning the high halo concentration 
c w 10 (concurring with the gravitational lensing analysis by 
Broadhurst et al. 2008; Lapi & Cavaliere 2009b) and involv- 
ing the flat slope a « 0.8 for the outer entropy profile. But 
we go beyond, and show in terms Eqs. (3) and (4) why these 
values deviate from the standard ones c » 4 and a k, 1.1, 
as spelled out above. In the same vein, CLFF09 find that a 
steeper density slope g = 3(a + bR)/5 w 2.4 is to apply in the 
outskirts. The high concentration of A1689 implies this to be 
an old structure with the bulk region dating back to a transi- 
tion epoch as early as « 1.5. In fact, the feature common to 
CC clusters like A2199, A2597, and A1689 is constituted by 
low values of a < 1 and high values c> 4 (see Table 1), that 
follow from their being generally old structures with shallow 
outer potential wells. 

At the other extreme, the NCC clusters appear to be dy- 
namically young structures from our determination of DM 



concentration and slope in the outer entropy profile. For ex- 
ample, in A1656 (Coma Cluster) our SM fit requires a high 
value a Pi 1.3, and relatedly (see Eqs. 3 and 4) a low concen- 
tration c « 4 and young age z, < 0.5. Toward the center, the 
NCC clusters are marked by a rising or flat temperature pro- 
file and by a generally flat brightness distribution. This occurs 
for central levels of kc exceeding some 50 keV cm^, and also 
points toward a thermodynamicaUy young age for the ICP. In 
fact, frequent and intense merger/AGN activity is expected in 
these clusters observed in the aftermath of their fast initial col- 
lapse, with considerable residual occurrence of mergers and 
AGN outbursts that lead to large central injections of energy 
and entropy. 

Such features are exhibited, in a sequence of increasing 
complexity, by A1656, A2256 and A644. Here, the SM ehc- 
its not only a high level, but also a pattern for the entropy 
deposited in the form of a floor extended out to rj. This we 
interpret in terms of the stallation radius attained by a power- 
ful, outbound blast either triggered by a major head-on merger 
(cf. simulations by Schindler et al. 2002, Vazza et al. 2009) 
or driven by a violent AGN outburst (see Forman et al. 2005; 
Cavaliere & Lapi 2006; McNamara & Nulsen 2007; Puch- 
wein et al. 2008), before being degraded into adiabatic sound 
waves of the kind caught in action by Fabian et al. (2006) in 
the Perseus cluster. To reach r/ « 250 kpc it takes a rather 
extreme merger deUvering about lO*"* erg and triggering a Se- 
dov blastwave that expands as Rs oc f'/^ f^/^ with decreasing 
Mach number; alternatively, it takes an AGN outburst of about 
10*^ erg continuously driving a blastwave to expand at con- 
stant Mach number with Rs oc E^l^t, see Lapi et al. (2005) 
and Cavaliere et al. (2006) for details. As discussed in § 4.4 
and 4.5, in the NCC clusters A1656 and A2256 analyzed here 
such values of r/ are accompanied by evidence of ongoing 
mergers, hallmarks of a recent cluster formation. 

This interpretation relates ry to the dating of the merger re- 
sponsible for the energy/entropy input; the good performance 
of the SM implies such a time to be intermediate between the 
blast transit time rf/Mi's^ 10"' Gyr (see CavaUere & Lapi 
2006), and the time 0.3(/tc/15 keV cm^)!-^ w 1 Gyr needed 
by radiative cooling to erode an entropy floor of about 50 keV 
cm^. 

Such a timing also guarantees that an accurate description 
of the ICP thermodynamic state for both CC and NCC clus- 
ters is provided by the SM based on the hydrostatic equilib- 
rium expressed by Eq. (1). To complete the issue, note that 
not only the equilibrium of the ICP is somewhat faster to at- 
tain than the DM's (see Ricker & Sarazin 2001; Lapi et al. 
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TABLE 2 

Fit parameters of Eqs. (A1) and (A2) 



a 


.« 


11 


1 


1.25 


0.750 


0.389 


21.465 


1.26 


0.756 


0.399 


14.374 


1.27 


0.762 


0.411 


10.614 


1.28 


0.768 


0.423 


8.316 


1.29 


0.774 


0.436 


6.751 



2005), but also that circularized data (integrated over annuli, 
see Snowden et al. 2008) tend to effectively smooth out local, 
limited deviations from spherical hydrostatics and to better 
agree with equilibrium. 

How does the SM face the challenge of complexity posed 
by substructures as observed in A2256 and A644? Interest- 
ingly, we still obtain good fits if we extrapolate the SM out to, 
or perhaps beyond its literal limits, toward conditions where 
T{r) varies on the scale r/, or differs around two locations 
(cf. Figs. 8-11); these conditions highlight the capabilities of 
the SM as a mere fitting tool. By the same token, the SM 
provides sharp snapshots of physical conditions even when 
these are spatially complex (as for A644), and strongly sug- 
gests that they may be traced back to two merger outcomes: 
a 'hot spot' imprinted in A2256 to partially erase a previous 
cool state and to yield a nascent NCC; or a 'cold drop' im- 
ported into the hot medium of A644 that will offset cooling. 
Focusing on the cold component, such cases may be termed 
as RCCs for Remnant of Cool Cores. 

In closing, we revert to Table 1 that provides an overall 
view of the SM parameters for the present cluster set, to 
stress the following points. First, columns 2, 3 and 4 visu- 
alize the sharp correlations among cluster classes (CC/NCC), 
basic properties of the DM halo (high/low concentration im- 
plying old/young dynamical age) and thermodynamic state of 
the ICP (low/high central entropy level, narrow/wide entropy 
floor, and steep/flat entropy slope). Second, columns 4 and 5 
visualize how from the density and temperature distributions 
over the full data range we obtain entropy profiles of com- 
parable quality to Zhang et al. (2008) and Cavagnolo et al. 
(2009). In particular, compared to the latter authors we find 
similar values of kc and somewhat lower values of a (with the 
upper bound on a discussed in § 2); we note that our values 
are derived by directly fitting the primary X-ray observables 
Sx{w) and T(w) with the SM, rather than deprojected profiles. 



Concerning virial radii, our values in column 7 systematically 
agree with literature evaluations; the agreements of the related 
virial masses are commented upon in § 4.3, 4.4, 4.5. Finally, 
column 8 visualizes the boundary values kBTu useful to intro- 
duce a discussion of the actual shock strengths related to the 
development of the outer halo and to preheating levels of the 
(Warm-Hot) IGM (see § 3). 

To sum up, our analysis of several NCC and CC galaxy 
clusters has shown how effective is the tool constituted by our 
Supermodel. On using a simple formalism and a fast algo- 
rithm, this leads to accurately fit in terms of a few physical 
parameters the many data points concerning diverse temper- 
ature and brightness shapes to constitute the library indexed 
by Table 1. Whence we extract sharp information concern- 
ing the evolutionary stage and the thermal history of clusters; 
in the inner regions these include the level (kc), pattern (rj) 
and timing of the entropy injections into the ICP, and in the 
outskirts the slope (a) of the entropy deposited by accretion, 
simply related to the halo concentration (c) and age (z,). 

Finally, the Supermodel offers predictions as for the den- 
sity and temperature profiles in the cluster outskirts out to 
the boundary r « 7? facing the IGM; for example, within a 
surrounding supercluster enhanced IGM preheating may be 
expected to weaken the accretion shocks and lower the tail 
of the brightness and temperature profiles. Such outer re- 
gions challenge most current instruments but are coming of 
age with SUZAKU (see Bautz et al. 2009; George et al. 
2009; Reiprich 2009), and will constitute a main target for the 
next generation X-ray telescopes planned to study low sur- 
face brightness plasmas, such as WFXT (see Giacconi et al. 
2009; also htt p : //wf xt . pha . jhu . edu/| l. These will 
open up the way to use galaxy clusters as probes of the sur- 
rounding Warm-Hot IGM, an interesting perspective that we 
will pursue in a following paper. 

Note added in proof: The evidence concerning rem- 
nants of cool cores in NCC clusters is also discussed from 
an observational point of view by M. Rossetti & S. Molendi 
(2009, A&A, submitted). 
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APPENDIX 

ANALYTIC FITS TO THE a-PROFILES 

Here we provide analytic fits to the radial runs of both the density and the circular velocity in the DM halos following the 
(isotropic) a-profiles. 

It is convenient to express these in terms of the popular and flexible parametric expression for generic density profiles 



r" 



l+wc" 
l+wicr)" 



(Al) 



■■r/R 



introduced by Zhao (1996) and constituting an extension of the empirical NEW (Navarro et al. 1997) formula (recall that r - 
is the adimensional radial coordinate, while c = R/r-2 is the concentration parameter). 

We cure the latter's unphysical features: diverging total mass and steep central cusp corresponding to angled potential, by 
deriving the parameters in Eq. (Al) directly from the Jeans equation and its derivatives, as shown by Lapi & Cavaliere (2009a, 
see their Appendix A). In fact, by differentiating Eq. (Al) it is easy to see that the latter approach implies w = -{2- s) /{2- s- qu). 
In addition, we find 

6 Kent 2 (70 -7a)^ 



5 7o-7a 3 



3Kcrit/5-2(7o-7a)2/3 ' 
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here 7a = 3 a/5, 70 = 6 - 3 a, 7;, = 3 ( 1 + q;)/2, and Kci-it(a) — 8.23 -AAA a. In Table 2 we report the parameter determinations for 
several values of a. 

Using Eq. (Al), the function v^(r) to be inserted in Eq. (2) of the main text may be explicitly expressed as 

-2.-. -2-. 2j^i [(3 - , g , 1 + (3 - ; -w (cr)"] 
" 2Fx[{'i-s)/u,q,\ + {'i-s)lu--wc"] 

in terms of the hypergeometric function 2Fu for a plot, see Fig. 1 in CLFF09. 

The resulting fits to the density and circular velocity profiles of the a-profiles hold to better than 15% in the whole range 
10"^ < cr < 10. In particular, for a = 35/27 = 1.296 we recover the exact solution of the Jeans equation found by Dehnen & 
McLaughlin (2005). 
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